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Abstract 

We present the full next-to-leading order corrected differential distribu- 
tions d'^a/dpT/dy, da/dpT and da/dy for the semi-inclusive process p + p —>■ 
H +' X'. Here X denotes the inclusive hadronic state and px and y are 
the transverse momentum and rapidity of the Higgs-boson H respectively. 
All QCD partonic subprocesses have been included. The computation is 
carried out in the limit that the top-quark mass — oo which is a very 
good approximation as long as mH,PT < 200 GeV. Our calculations reveal 
that the dominant subprocess is given hj g + g H +' X' but the reaction 
g + q{q) H +' X' is not negligible. Another feature is that the i^-factor 
representing the ratio between the next-to-leading order and leading order 
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differential distributions is large. It varies from 1.4 to 1.7 depending on the 
kinematic region and clioice of parton densities. We sliow tliat a reliable 
determination of the differential cross sections requires good knowledge of 
the gluon density in the region where x < 10~^. Further we study whether 
the differential distributions are dominated at large transverse momentum by 
soft-plus- virtual gluon contributions. This is of interest for the rcsummation 
of large corrections which occur near the boundary of phase space. We also 
compare our results with those previously reported in the literature. 

PACS numbers: 13.85.-t, 14.80.Bm. 
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1 Introduction 



The Higgs boson, which is the corner stone of the standard model, is the 
only particle which has not been discovered yet. Its discovery or its ab- 
sence will shed light on the mechanism how particles acquire mass as well 
as answer questions about supersymmetric extensions of the standard model 
or about compositeness of the existing particles and the Higgs boson. The 
LEP experiments [1] give a lower mass hmit of about uih ~ 114 GeV/c^ 
and fits to the data using precision calculations in the electro-weak sector of 
the standard model indicate an upper limit rriH < 200 GeV/c^ with 95 % 
confidence level. After the end of the LEP program the search for the Higgs 
will be continued at hadron colliders in particular at the TEVATRON and 
the LHC. If the Higgs mass is in the above range the principal production 
mechanisms at hadron colliders are gluon-gluon fusion g + g ^ H +' X' or 
W^+W^~-fusion appearing in the reaction q + q^q + q + H+' X'. In these 
processes 'X' denotes an inclusive hadronic state. For mn < 135 GeV/c^ the 
process q + q^V^H + V+'X' with V — '-f,Z,W has also to be taken 
into consideration. The Higgs boson will be observed via its decay products 
among which the channels H ^ b+b and H ^ ^ + '~f are the most prominant 
ones although if — > r"*" + r~ should also be considered. However the large 
backgrounds make the observation of these decays very difficult and it will 
take a lot of experimental and theoretical effort to detect the Higgs boson 
provided it is there. 

In this paper we concentrate on Higgs production channels where the 
lowest order reaction proceeds via the gluon-gluon fusion mechanism. In the 
standard model the Higgs boson couples to the gluons via heavy quark loops 
among which the top-quark loop is the most prominent since the coupling 
of the Higgs to a fermion loop is proportional to the mass of the fermion 
(for a review see [2]). In lowest order (LO) the gluon-gluon fusion process 
g + g ^ Hj represented by the top-quark triangle graph, was computed in [3]. 
The next-to-leading (NLO) processes given by gluon bremsstrahlung g + g —>■ 
g + H and g + q{q) — > q{q) + H were presented in [4], [5] and [6] from which 
one can derive the transverse momentum {pt) and rapidity {y) distributions 
of the Higgs boson. The total integrated cross section, which also involves 
the computation of the QCD corrections to the top-quark loop, has been 
calculated in [7]. This calculation is rather cumbersome since it involves the 
computation of two-loop triangular graphs with massive quarks. Recently 
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also the next-to-next-to-leading (NNLO) processes involving all two-to-three 
parton processes have been computed in [8] using the helicity method which 
means that the matrix elements are presented in four dimensions. Prom the 
experience gained from the NLO corrections it is clear that it will be very 
difficult to obtain the NNLO one-particle inclusive distributions from these 
calculations let alone the total cross sections. 

Fortunately one can simplify the calculations if one takes the large top- 
quark mass limit ^ oo. In this case the Feynman graphs are obtained 
from an effective Lagrangian describing the direct coupling of the Higgs boson 
to the gluons. The LO and NLO contributions to the total cross section 
in this approximation were computed in [9] and they found that the error 
introduced by taking the nit ^ oo limit is less than about 5% provided 
TTiH < 2 m^. A similar investigation was done for the differential distributions 
of the processes g + g ^ g + H and g + q{q) — > q{q) + H in [10]. Here the 
approximation is valid as long as rriH and are smaller than m^. This is 
corroborated by the recent calculations in [8] which show that for itlh = 
120 GeV/c^ the approximation is valid for jets with transverse momenta 
smaller than 200 GeV/c. The NNLO matrix elements using the effective 
Lagrangian were computed in [11], [12] albeit in four dimensions. The one- 
loop corrections to the two-to-two parton subproccsscs were presented in [13], 
where the computation of the loop integrals was performed in n-dimensions 
but the matrix element was still presented in four dimensions. The results 
in [11], [12] and [13] were used to compute the transverse momentum and 
rapidity distributions of the Higgs boson up to NLO [14]. The effective 
Lagrangian method was also applied to obtain the NNLO total cross section 
by the calculation of the two-loop corrections to the Higgs-gluon-gluon vertex 
in [15], the soft-plus- virtual gluon corrections in [16] and the computation of 
the two to three body processes in [17]. 

In this paper we present the full NLO computation of the double differen- 
tial distributions dPa/dpT/dy for Higgs boson production in hadron-hadron 
collisions using the gluon-gluon fusion mechanism in the ^ oo approxima- 
tion. Here we have included all partonic subprocesses. A similar calculation 
has been performed in [14] but our approach differs from it in various as- 
pects. First our calculation is purely analytical and follows the calculation 
carried out for the Drell-Yan process describing vector boson production in 
hadron-hadron collisions (see [18]). The approach in [14] was mainly nu- 
merical and based on the methods explained in [19]. Moreover it used the 
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two-to-three particle matrix elements in [11], [12] and the two-to-two particle 
matrix elements including virtual corrections in [13] which were all presented 
in four dimensions. In our calculation the matrix elements as well as the 
loop integrals and phase space integrals are computed in n dimensions. The 
advantage of the analytical approach is that one gets more insight into the 
structure of the radiative corrections. This is particularly important for the 
large corrections, due to soft gluon radiation and coUinear fermion pair pro- 
duction, which arise near the boundary of phase space, where the pr of the 
Higgs boson gets large. Resummation of this type of corrections has been 
carried out for the total cross section in [20]. Resummation of small pt con- 
tributions due to the Sudakov effect has been done in [21]. In view of the 
experimental problems to observe the Higgs boson, a recalculation of all the 
NLO corrections is necessary to be sure that the theoretical predictions are 
correct. Wc find that the contribution due to the (anti-) quark gluon subpro- 
cess is substantial in particular at large transverse momentum and that the 
important region of x in the gluon density is a; < 10~^. We also investigate 
the region of applicability of the soft-plus- virtual (S-l-V) approximation for 
the calculation of the differential cross sections. Finally we mention that an- 
other paper has just appeared on the NLO corrections to the g + g— > H + g 
channel, using the helicity framework [22]. 

Our paper will be organized as follows. In section 2 we give an outline of 
the kinematics and present the Born contributions. In section 3 we present 
the virtual contributions. In section 4 the gluon bremsstrahlung corrections 
to the Born reactions are computed and the soft gluon cross sections are 
explicitly shown. In section 5 we show how the mass factorization is car- 
ried out. In section 6 we give differential cross sections for proton-proton 
collisions at the LHC and make comparisons with results obtained earlier in 
the literature. Some particular expressions are given in appendices A and B. 
Other formulae, in particular those for the two-to-three body cross sections, 
are too long to be published. They are available upon request as files written 
in the the algebraic manipulation program FORM [23]. 
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2 Lowest order contributions to Higgs pro- 
duction 

In the large top-quark mass hmit the Feynman rules (see e.g. [12]) can be 
derived from the following effective Lagrangian density 

C,ff = G^{x)0{x) with 0{x) = -\gI^{x)G^^^'^{x) (2.1) 

where ^{x) represents the Higgs field and G is an effective couphng constant 
given by 

G^ = 4V2 (^)' G,t'F\t)C' ■ (2.2) 

In the expression above m and mt denote the masses of the Higgs boson and 
the top quark respectively. The running coupling is given by o^s{^J^) where 
IXr denotes the renormalization scale and G^? is the Fermi constant. Further 
C is the coefficient function which originates from the QCD corrections to 
the top-quark triangle graph describing the process H ^ g + g ui the limit 
mt — > oo. On the Born level the width of this decay process is given by 

"^(H^g + g)- o^M) F\r) , (2.3) 

where N denotes the number of colours. Further the function F{t) occurring 
in Eqs. (2.2), (2.3) is defined by 

F(r) = 1 + (1 - r) /(r) , r = ^, 
/(r) = arcsin^ , for r > 1 , 

/(r) = -i(lnl^^ + .z)' for r<l. (2.4) 
In the large mt-limit F{t) behaves as 

limF(r) = ;^. (2.5) 
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The coupling G in Eq. (2.2) is presented for general mt on the Born level only 
whereas C is computed in higher order for mj — > cxo. In order to keep some 
part of the top quark mass dependence we take for G the expression in Eq. 
(2.2). This is an approximation because the gluons which couple to the Higgs 
boson via the top-quark loop in the partonic processes describing Higgs- 
production arc often virtual. The virtual-gluon momentum dependence is 
neither described by F{t) in Eq. (2.4) nor by C. For on-mass-shell gluons 
the latter quantity has been computed in the large top-quark mass limit up 
to order ag in [7], [9], [11] and up to in [20], [24]. In second order it reads 



C as{i4), 



j4_ 




18 rrit 



/ 67 8 u2 ^ 



(2.6) 



Here Uf denotes the number of light flavours and af^ is presented in a five 
flavour number scheme. 

In this paper we study the semi-inclusive reaction with one Higgs-boson 
H in the final state. It will be denoted by 



(2.7) 



where Hi and H2 denote the incoming hadrons and X represents an inclusive 
hadronic final state. In lowest order the partonic reactions contributing to 
Eq. (2.7) are denoted by 

a(Pi) + b{p2) c(-p3) + H{-p^) , Pi + P2 + P3 + P5 = , 

a,b,c = q,q,g, (2.8) 
and the partonic cross sections in n-dimensions are given by 



a 



(1) 

ab^c H 



^ ^ 2s J (27r)"-i J (27r)"-i ' 



(2.9) 



Here M^^^ denotes the amplitude of the process and the strong coupling con- 
stant is given by g with g"^ = Anag- The scale 11 originates from the fact that 
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the coupling constant acquires a dimension in n dimensions. The quantity 
Kab represents tlie spin and colour average over the initial states including 
the statistical factor 1/m! if one integrates over m identical particles in the fi- 
nal state. Further we have assumed that the Higgs boson is mainly produced 
on-mass-shell. However one can also use the narrow width approximation. 
This can be achieved by replacing 

^ipl-rn') - (2-10) 

where F is the total width of the Higgs boson which is dominated by the 
decay H ^ g + g. Choosing the CM. frame of the incoming partons we have 
the following parametrization 



Pi 



Jv^i2(l,0,---,0,1), 



P2 = ^V^12(l,0,---,0,-l), 



-P3 = -7—7^ — (l,---,-sm^i,-cos^i), 
2v -r 12 



-P5 = (P12 + • • • , (P12 - m^) sin 61, {P^ - m^) cos ^1) , 
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(2.11) 



with 

Pij = {Pi + Pjf , -P12 = s , Pib^t, P25^u, s + t + u^m^ , 
t — u 

cos^i = -. (2.12) 

s — 

From Eq. (2.9) we infer that 

dtdu r(l-F£/2) 47r \lJp s J 

x\Mj^l^H\\ with n = 4 + e. (2.13) 
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The spherical factor is defined by 

'5e = expQ(7E-ln47r)^ . 

On the Born level we have the following subprocesses 

9 + 9^g + H, 



K. 
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\m: 



(1) 



qq^g H\ 



(1) 



qg^q H\ 



stu 



+2 £171 stu 



4(l + £/2)2(Ar2_l)2 ' 

q + q ^ g + H , 
1 



Ca Cf 
1 1 

q{q)+9 ^ q{q) + H, 
1 

CaCf — 



1 + -] {r + u'} +etu 



Kqg = 



1 



u 



1 



est 



4(1 + e/2) iV (iV2 - 1) ' 



with 



Ca^N, Cf^ 



The n-dimensional matrix element |M^^^^ hV '^^ proportional to n — 2 pro- 
vided m = so that it vanishes in two dimensions. The same feature also 



2N 
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appears in the lowest order matrix element for the decay H g + g de- 
rived from the effective Lagrangian in Eq. (2.1). Notice that the factor 
n — 2 — 2{1 + e/2) also shows up in the spin average quantities Kab- 
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3 One-loop corrections to the 2 2-body 
reactions 

The one-loop corrections to the gluon-gluon subprocess in Eq. (2.15) entails 
the computation of forty-two graphs. Fourteen graphs lead to independent 
expressions whereas the remaining ones can be obtained via crossing from 
the s- to t-channel or from the s- to -u-channel. The Feynman integrals show 
ultraviolet, infrared and collinear singularities which will be regularized using 
n-dimensional regularization. Hence the matrix element and the Feynman 
integrals have to be computed in n dimensions for which we used the alge- 
braic manipulation program FORM (version 3.0) [23]. Since the Feynman 
integrals also contain the integration momentum in the numerator we have to 
apply tensorial reduction. Here we followed the procedure in [25] which was 
extended to n-dimensions in [26]. The scalar integrals for the two-, three- 
and four-point functions can be found in [27] . The one- loop correction to the 
gluon-gluon differential distribution (2.15) can be written as 
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(1) I 

qq-^g H\ 



(3.3) 



\MB^^^^H?-CACF{-t-u). (3.4) 

The square of the matrix element for the virtual corrections to the quark- 
gluon subprocess in Eq. (2.21) can be obtained from the one calculated for 
the quark-anti-quark subprocess in Eq. (2.18) via crossing from the s-channel 
to the li-channel. This yields 




(g — ^2 \ / t \ t U 
+Li2 — -21n — In — 
s J J jjL'^ iJ,^ 
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^, —u, s ^, n —u 1 r, t — m 
+21n — In— -21n2 — - -In^ — 



+;in^ 



m —t 



-t 



—u 5 



+C(2) - 
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IM, 



(1) 



qg-*q H\ 



+ {Ca-Cj.){^-}\MB[ 



,(1) I 

'qg^q Hi 



(3.5) 



|MsJJL,/^|' = CaC^(« + ^)- (3.6) 

Notice that the MB^^^ in Eqs. (3.2), (3.4) and (3.6) are presented in four 
dimensions and higher order terms in e are ignored. This is sufficient because 
the factors with which they are multiphed are finite in the hmit £ — > 0. Our 
results agree with those presented in [13]. Notice that in the latter paper 

the Born amplitudes M^l^^ ^ are presented in four dimensions because they 
are constructed using the helicity method. We have computed them in n 
dimensions. Apparently this does not affect the factors that multiply the 
Born matrix elements and which contain all the pole terms in 
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4 Gluon bremsstrahlung and other 
2 3-body processes 

In this section we give an outline of the computation of the two-to-three 
parton processes which show up in NLO. The calculation proceeds in an 
analogous way as in the case for heavy flavour production presented in [29]. 
The processes under consideration will be denoted by 

5 

a{pi) + b{p2) c{-ps) + d{-p4) + H{-p5) , ^p. = , 

i=l 

a,b,c,d = q,q,g, P • = , i = l-4, pl = m^. (4.1) 

The cross section corresponding to the reaction above can be expressed in n 
dimensions as follows 

.(2) _ ^ J_^2 4 [ ^>3 [ dy^ r d^r, 

(^ab^cd H 2s ^ i (27r)-i J (27r)«-i J (27r)"-i 

^'-{pI) 5+{pI) S^ipl - m^)(27r)"5(")(pi +p,+p,+p,+p,) 

x\M^lcdH{0i,02)\\ (4.2) 

where Kab is defined below Eq. (2.9). Choosing the CM. frame of the out- 
going partons 3 and 4 we get the following parametrization in n dimensions 

pi=u;i(l,0,---, 0,0,1), 

P2 = (<^2,0, • • • ,0, |P5| Sin^, |p5| COSV' - UJi) , 

—ps = a;4 (1, 0, • • • , sin 9i sin $2, sin 9i cos 02, cos ^i) , 

— P4 = a;4 (1, 0, • • • , — sin 9i sin 6*2, — sin 9i cos 6*2, — cos ^i) , 

-P5 = (^^5, 0, • • • , 0, IpsI sin V', IpsI cos^') , (4.3) 
where the energies uji arc given by 
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Pi5 + P25 I Pi5 - m"^ + 2ui LJ5 
UJ5 — 1= — , cosw = —r, . (4.4 

The invariants are denoted by 

Pij = {pi + PjY , -P12 = s , Pi5^t, P25 = ^i , 

P34 = S4 = s + t + M - . (4.5) 

Prom I cos'^l < 1 and S4 = P34 > we derive the boundary conditions 

2 

<m^-u<s + t, <m^ -t<s. (4.6) 

—t 

In the subsequent part of this paper it is more convenient to choose S4 as 
integration variable instead of u. In this case the integration boundaries are 

tifn^ — s — t] 

< S4 < ^ ^ , <m^ -t<s. (4.7) 

— t 

Prom the above kinematics the differential cross section in Eq. (4.2) becomes 
equal to 

dtdu 2 "'r(l + £) I 47r M s i W) 



x|Mi^L,,^|^ (4.8) 



where |Af^j,i^^.^| //P is the second order matrix clement integrated over the 
polar angle 61 and the azimuthal angle 62- It is given by 



iM^lca h\' = r de, (sin 0,y+^ r d02 (sin e^Y \M^l^ ^{0,, 9, 
Jo Jo 

^ j dn^_,\M^^l,, ^{OM? . (4.9) 

As in the calculation of the lowest order matrix elements in Eqs. (2.16), 
(2.19) (2.22) the next-to-leading order expressions |M(2)(^i,^2)P have to be 
computed in n dimensions because we use n-dimensional regularization for 
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the collinear and infrared singularities which arise from the integration over 
the momenta of the final state partons. We use the algebraic manipulation 
program FORM (version 3.0) as in the earlier computation of the virtual 
corrections. Further we choose the Feynman gauge for the gluon propagators 
and the sum over the physical polarizations of the external gluons is given 
by 

J2 e^{p,a)e''{p,a)^P^'''{l,p), with P'''' ^ P^" I, ^ , f^O, 

a=L,R 

(4.10) 

with 

I ■ p 

where / is an arbitrary lightlike vector. There are three types of matrix 
elements depending on the number of gluons and (anti-) quarks appearing 
in the initial and final state. They will be denoted by ggggH, ggqqH and 
qqqqH . Notice that the quark q can also represent an anti-quark q depending 
on the type of process. The Feyman graphs can be found in [12]. The matrix 
element for ggggH involves the computation of twenty-six graphs whereas 
eight graphs contribute to ggqqH. In the case of qqqqH one has to distinguish 
between identical and non-identical quarks. For identical quarks we have two 
Feynman diagrams whereas the non-identical quarks are represented by one 
graph only. In the case of ggggH the amplitude is denoted by M'^'^^^'^'^p 
and the matrix element squared is obtained by contraction over the Lorentz 
indices as follows 

\M^'\e-,,e,)\^ = P'^'^P^^P'^^P^-M^XpM^X ■ (4-12) 

We checked that the dependence on I -pi disappears for arbitrary /. However 
this was only possible by writing the expression in Eq. (4.12) over a common 
denominator and expressing the result into five independent kinematical in- 
variants. It turns out that Eq. (4.12) has the same property as the square 
of the Born matrix clement in Eq. (2.16) namely that it is proportional to 
n — 2 provided m = 0. In the case of ggqqH the amplitude is given by M^J . 
Here we checked that the dependence on / • of the matrix element squared 
given by 

|M(2)(^i,^2)r = P'"' P"^ Mj^J M^^^ , (4.13) 
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already disappears after a simple partial fractioning. We made a compar- 
ison with the results in [12] which are presented in four dimensions only 
because they are computed using the helicity method. After corrections for 
some misprints ^ the expressions obtained from Eqs. (4.12) and (4.13) are 
in agreement with those presented in Appendix A of [12]. The procedure to 
bring the matrix clement over a common denominator, as discussed below 
Eq. (4.12), is not suitable for partial fractioning (see below) because it leads 
to high powers in the angular dependent kinematical variables, which show 
up in the numerator and denominator. To avoid this complication we proceed 
in a different way. Since the l.pi terms are linear independent for i = 1 — 4 
they have to cancel among themselves. This we checked by deleting them 
and comparing the truncated matrix element squared with the one obtained 
by bringing all terms over a common denominator, which is manifestly in- 
dependent on I. Pi. The difference between these two expressions turned out 
to be zero. Moreover the total power of the angular dependent kinemati- 
cal invariants appearing in the numerator and denominator of the truncated 
matrix element does not exceed four so that it becomes amenable for partial 
fractioning. After this procedure the integration over the angles is performed 
in the following way. The expression \M^'^\9i, ^2)^ depends on the variables 
Pij defined in Eq. (4.5) among which five are linearly independent. Because 
of the parametrization in Eq. (4.3) the variables -P12, -P15, -P25, -P34 are inde- 
pendent of the angles 9i, 62 whereas the remaining ones P13, P23, -P14, -P24, -P35 
and P45 are angular dependent. One can distinguish two types of integrals 

J(i^7^PJ) = j d^n-iPr'P^i,, ^,J = l-4 m = 3,4, 

with -4<A;<2 0</<2, (4.14) 

with -3<k<2 -3<l<2. (4.15) 
The basic integrals of the first class are given by 

JjPrs^PAB) - 2 ^ (P25 - Pu)-' Ik,i{A B, C) , 

^We found additional misprints in [12]. The 534 in the denominators of the two terms 

in (A. 12) should read ^24. In (A. 16) one has instead of ni3 = —5*24 «i2(2 ^ 4) and 
n23 = —523 f^i3(3 ^ 4) the expressions ni3 = —5*24 ni2(l ^ 4) and 7123 = ni3(3 ^ 4). 
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J{Piz\Pz^) = 2^{P2,-PM)-'h,i{A-B,-C), (4.16) 

with 

Ik,i{A^B,C) = J dnr,-i{l - cosOi)-'' {A + B cosOi + C sinOi cos^s)"'- 

(4.17) 

The basic integrals of the second class are represented by 

J{Plt,P2A) = 2^+'(P25-P34)-'(Pi5-P34)-'/m(x), 

J{Pn,P2j) = 2'=+^(P25-P34)-'(Pi5-P34)-^/m(^-x)> (4-18) 

with 

h,i{x) = y^(in„_i(l — cos6'i)~'^ (1 — 0086*1 cosx — sin6'i cos6'2sinx)~' 
a-fc-^ r(n/2 -1-k) r(n/2 -l-l)T{n- 3) 



= 7r2^ 



xF,,2[k,l,^ -l,cos'^j , (4.19) 

where -^1,2(0, b, c; x) denotes the hypergeomctric function which can be found 
in [30]. The other integrals can be derived by interchanging the inclusive 
momenta pa and or interchanging the parametrization for pi and p2- The 
angular integrals above are presented in Appendix C of [29] for the case 
—2 < k <2. Because the square of the matrix element for Higgs production 
contains two extra powers of Pij in the numerator the previous computation 
has to be extended to cover the cases k = —3, —4. These extra powers can be 
attributed to the momentum dependence in the effective Higgs-gluon-gluon 
coupling described by the Lagrangian in Eq. (2.1). The colour decomposi- 
tions of the integrated expressions read as follows 

g + g^g + 9 + H, (4.20) 



g + g^Qi + qi + H, (4.22) 
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(2) 



= UfCACp 



(2) 



gq^gg H\ 



q + q^g + g + H, 
Ca Cp 



(4.23) 
(4.24) 



(2) 



(2) 



gg-*qg H\ 



q + g ^ q + g + H , 
Ca Cf 



CA\Mill,,H\\+CF\M, 



(2) 12 
qg-*qg Hi f 



, (4.25) 
(4.26) 
. (4.27) 



The colour and spin average factors Kgg, Kgg and Kgg are given in Eqs. (2.17), 
(2.20) and (2.23) respectively. If two identical particles in the processes 
above or below appear in an inclusive final state a factor 1/2 is implicitly 
understood. In NLO we encounter some new subprocesses. The first one 
is given by the reaction in Eq. (4.22) where a sum over all light flavours 
qi with i = 1- ■ - Hf is understood. The next one is quark-quark scattering 
(non-identical and identical quarks) represented by 



qi + q2 ^ qi + q2 + H , qii^qi-, 



\M, 



(2) 



Ca Cf Wq^q^^q^qi H?' A' 

q + q^q + q + H, 



|M. 



(2) 



CaCf\M;,±_^^^^^^\\ + Cf\M^'^ 



gg—*ggH\ ^a^p i"'^gig2— >?ig2 -hi a 

which have the colour and spin average factors 

1 1 



qq^qq H\ f 



K. 



K = 

4 7V2- 



The second subprocess is quark-anti-quark scattering 

qi + q2^ qi + q2 + H , qi q2 , 



\<i-^q.q.H\' = CaCf\M^'^ 



?192— >?iq2 H\ A^ 

q + q^qi + qi + H, qi 7^ q , 



(4.28) 

(4.29) 
(4.30) 

, (4.31) 
(4.32) 

(4.33) 

(4.34) 
(4.35) 
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q + q^ q + q + H , (4.37) 



l^if-^qq i? 1^ - CaCf i l-^Sa^gigj H?A + \^qf-^qaqa h\'^ A 



+Cf\M^II,,h\'f^ (4.38) 
and the colour and spin average factor reads 

K,,,, = K,, = . (4.39) 

Further we have summed over aU quark flavours in the final state provided 
they arc not cqiial to the flavour of the quark q in the initial state (sec Eq. 
(4.35)). Notice that the colour decomposition given above also holds before 
the squares of the matrix elements are integrated over the angles. Some of the 
expressions in Eqs. (4.20)- (4.38) are singular in the limit S4 0, where S4 is 
defined in Eq. (4.7), which is due to soft gluon radiation or soft (collinear) 
fermion pair production. This will lead to infrared singularities when the 
differential cross section in Eq. (4.8) is convoluted with the parton densities 
even after mass factorization is carried out. The convolution involves an 
integration over S4 so that one has to split up the partonic cross sections into 
hard and soft gluonic parts as follows 

j2 -(2) .2 -HARD j2 ^SOFT 

dtdu dtdu dtdu 

(4.40) 

where A is chosen in such a way that it satisfies the inequalities 

A<s, A^-t, A^-u. (4.41) 
The soft gluon contribution to the differential cross section is defined by 

j2 -soft 1 02 / /',,2U 2 / + „ \ £/2 

2 " ^ab-*cd H _ ^ '^e .^2 / ^M/^ 

dtdu 2 "'r(l + £) 1^ 47r 

5(s + t + u-m^) j^ds^ (^) \M'^^^jj 




2 



(4.42) 
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Only the singular part of iM^^i^^^ ji^p, which behaves as I/S4, contributes to 
the above integral whereas the non-singular terms vanish in the limit S4 ^ 0. 
This singular part is called \M^^J^J^ and behaves in two different ways 

I. 

WT!^ - I , (4.43) 

and 



II. 

WSI^-li^f- (4.44) 

The (Pa^^^'^ for the various subprocesses are presented in Appendix A. 
Adding the latter to the virtual contributions presented in section 3 leads 
to the soft-plus-virtual differential cross sections defined by 

J2-S+V J2-SOFT r/^n-VIRT 

at a u at a u at a u 



The infrared singularies cancel in this expression so that the double pole 
terms occuring in both the soft and the virtual contributions, all vanish. 
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5 Renormalization and mass factorization of 
the partonic differential cross sections 

The virtual cross sections in section 3 have two types of ultraviolet diver- 
gences. One is removed by renormalization of the operator 0{x) occuring in 
Eq. (2.1) and the other by renormalization of the strong coupling constant. 
The renormalization constant corresponding to the former is 



1 + 

e 



a.s 2 



Att e 



1 + Z^'^ + Z'n' + 



a. 



47r 



+ 



(2) 



O 



(5.1) 



up to order «^ [31], where /?o and (3i are the lowest order coefficients of the 
bet a- function given by 



47r 



A-kJ 



/5o — 
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10 



Pi = -CX-2nfCF--nfCA. (5.2) 



In lowest order Zq is the same as the renormalization constant for the strong 
coupling constant which reads 



a. 2 



+ 



(5.3) 



Since the operator insertion always appears twice in the differential cross 
sections the latter have to be multiplied by Zq. 

After renormalization one has still to perform mass factorization to re- 
move the remaining coUinear divergences. This is achieved by the formula 



d t d u 



iu Jo Xi Jo Xo 
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a t a u 

where an implicit summation over e, / is understood, with 

S — X1X2S, t — Xi{t — m^) + , u — X2{u — m^) + . (5.5) 

In Eq. (5.4) cPa represents the singular cross section which still contains 
the coUinear divergences indicated by £. These divergences are removed by 
the kernels Tab so that the finite cross sections d^a are left. Both therefore 
depend on the mass factorization scale /i. Since the mass factorization has 
to be carried out for the semi-inclusive cross section up to NLO the kernels 
Fab have to be corrected up to order ag only. They are denoted by 

Tabix^s) = SabS{l -x)+ r'^^x^e) . (5.6) 

In section 4 the cross sections were split up into hard gluon {H) (Eq. (4.40)) 
and soft-plus-virtual gluon parts (S'+K) (Eq. (4.45)). Therefore one proceeds 
with the kernels in the same way so that the mass factorization can be carried 
out for both parts separately 

ril^ = rr°'^^) + 5(i-x)rr'^^ (5.7) 

The various kernels are expressed into the splitting functions Pab, depending 
on the partons a,b, which represent the residues of the coUinear singularities 



11 


J- qq 




r(i) 

91 


_ r(i) 

91 




r(i) - 

19 


_ r(i) _ 

^ 19 





The splitting functions are given by 
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'1+ (1 -a;)^ 



X 



1 



1 11 

, +- -2 + x-x^ + — (5(l-a;) 
1 — X 12 



--nf6{l-x) . 



(5.9) 



The coUinearly finite cross sections can be derived from Eq. (5.4) by using 
the above kernels. In LO they arc equal to the Born results derived in section 
2. In NLO the finite order contributions can be written as 

9 + 9^9 + 9 + H, 



^2 " ^gg^gg H 
d t d u 



— s 



d^ a^^^ 
2 " "gg^gg H 

dtdu 



^ ^2 Q (^gg^g H _ .2 " ^gg^g H 

didu didu ' 



(5.10) 



g + g^q + q + H, 



d? a^'^ 
dtdu 



— s 



2 " ^gg-^qq H 

dtdu 



d^ a^^^ d^ a^^^ 

-2 rw ® 5^-^ ^ - 2 rw .^^-'^ ^ , 

dtdu dtdu 



(5.11) 



q + q^g + g + H, 



j2 " ^qq^gg H 
dtdu 



— s 



d^ a^^) 
2 " "qq^gg H 

dtdu 
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j2 (1) j2 (1) 

dtdu dtdu 



(5.12) 



q + 9 ^ q + 9 + H , 



d^ a^^^ r/2 ^(2) 

dtdu dtdu 



j2 (1) ^2 (1) 

yy ri -f- HI ^ ^ 



dtdu dtdu 

_r(i) ^ ^2 _ rw ^ ^^2 

dtdu dtdu 

(5.13) 

gi + g2^gi + g2 + -f^, q + q ^ q + q + H , 



r/2 rr^^) r/^ ^^^^ 

dtdu dtdu 



j2 (1) ^2 (1) 

_r(i) ^ ^2 _ rw ^ ,^2 ^ 

dtdu dtdu 

(5.14) 

gi + ?2^gi + ?2 + -H", gr + ^^g + ^+ i?, 



rf2 a^^) rf2 ^(2) 

„2 " "qq^qqH ^ ^2 " " qq^qq H 

dtdu dtdu 



j2 (1) j2 (1) 

_p(i) ^ 52 ^ ^gg^g-tf _ -pa) 5> 52 ^ ^gg^g h 
dtdu dtdu 



(5.15) 



Notice that the finite parts of these cross sections can become negative due 
to the above subtractions. The expressions for the finite hard gluon cross 
sections d"^ Cab^d h too long to be published. They exist as FORM 
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files and they are available on request. Here we can only show those parts 
which behave like I/S4 since they are of interest later on. In the soft gluon 
limit S4 — > they read 



^2 ^HARD 
lim 99^99 

s4^o a t d u 



N 



1 



47r J {m - 1)2 S4 



31n — 



-In 



tu 11 



\M. 



(1) 



gg^g H\ i 



(5.16) 



j2 ^HARD 

54_*o a t a u 



77 G 



2 / "^(/^^)\^ 



47r y (Ar2 - 1)2 



^2 ^HARD 

lim — gg-^gg ^ 
s4-^o a t a u 



xlM, 



(1) 



gg^g ^1 ' 

2N\ 2 



47r i 5^ 



C . - In 



£4 



(5.17) 



, tu 11 

+ In 

fi^s 12 



CF<41n^ -21n 
I 



(5.18) 



j2 ^HARD 

s4-*o at a u 



47r J Ar(7V2 - 1) S4 



.l„i!ii+cJln4-- 



IM^^L., ^1^(5.19) 



^2 ^HARD 
lim ^gg->gg-^ 
s4->o at du 



47r 



|M, 



(1) 



(5.20) 



Notice that the expressions above satisfy a supersymmetric relation. It turns 
out that the expressions for gg X' + if in Eqs. (5.16) plus (5.17), qq — 
X' + H in Eqs. (5.18) plus (5.20) and qg ^' X' + H in Eq. (5.19) become 
equal for a jV = 1 supersymmetry where Ca = Cp = Uf = N. 



28 



After renormalization of the the operator 0{x) and the couphng constant 
using Zo and in Eqs. (5.1) and (5.3) respectively the ultraviolet di- 
vergences vanish in the soft- plus- virtual cross section (P a^^^ in Eq. (4.45). 
The remaining collinear divergences are removed via mass factorization using 
T^^^ in Eqs. (5.7)-(5.9). The combined effect of these three operations leads 
to the following expressions 



dtdu 



— s 



i " '^gg^g H 

dtdu 



+ Z« + 2 Z^'^ 



o 



a 



(1) 

gg^g h 



dtdu 



(5.21) 



J2 S+V 
! Q ^qq^g H 

dtdu 



qq-^g H 



dtdu 



(1) 



d'' rr^^) 
dtdu 



(5.22) 



d"^ rTS+^ 
dtdu 



— s 



dtdu 



d- a^^ 



d' a^y 



dtdu 



gg 



dtdu 



_pS+V,(l) 2 " "qq^g H 

«^ dtdu ' 



(5.23) 



The results for the finite soft-plus-virtual differential cross sections read 
TT 5(s i u - m^) 



d'^ ^7^+^ 
.2" ^gg^g H 



1 



dtdu 



47r 



(iV2 _ i)i 



X 




+7V -ln2-+(-ln— --)ln — + L12 
\2 //^ ^ //^s 12^ /x^ 
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f u \ ^ ( s — m 

+Ll2 ^ + Ll2 



In ^ In 4^ - In In ^ 



1 , 9 — t 1 o — M 1 9 s 1, r, f t — m^' 

+- In^ — + - In^ — + - In^ — In^ 

2 2 ij? 2 2 \ t 



--In^ +-ln2 — +-ln2 — 



+ {N-nf){\}\MB^^ 



(1) 12 

gg-*g HI 



(5.24) 



dtdu 



47r / 



' A, A 1 ,9 
n/ I - in — + - m 



tu A ^. s-m^ 
1 In — + L12 



\^ 2 //^ ^ //^s 12' //^ V s 



-In' li^2 , 1,^2 5 11^ s 73 ~ 



--In' — + -ln^ — -— In— -5C(2) + 



/x^ 6 /x^ 



12 



f 2 In' 4 - 2 In 4^ In 4 + Li2 f 4 1 + 



IJ?S jJ? 



u 



In — In — - In — In — + In — + In — 
/x^ /x^ /x^ /x^ /x^ /x"^ 
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4ln^f^l4ln^ + 8C(2)-4))|Ae,„l^ 



dtdu 



7r5(s + t + u-m'^) G 



(5.25) 



47r J N{m - 1) 



X 



tu , A 



-M , s 



+Ca In^ — - In — In — + Li2 ^ - In — In — 
+ ln' — --InM +-lnM ^ 



2 2 



1^ 



+4C(2) + f 



— In — 
6 //2 



+Ci. I ^ In^ 4 - T In 4 + Li2 ( ^ ) + LI 



.2 1^' 4 z^' 



s — m 



, -t , S 1 , 2 ~^ 1 1 2 1,25 



In^ I : I + - InM — + - In 



t 



2 //^ 



The expressions multiplying the Born matrix elements \M^ll^^ ji^p satisfy 



(5.26) 
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a supersymmetry relation. If wc choose a, M = 1 supersymmetry so that 
the quarks are put in the adjoint representation i.e. Ca = Cp = Uf = 
N then all expressions are equal except for the rational constants which 
arc equal to 19/12 in Eq. (5.24), 15/12 in Eq. (5.25) and 17/12 in Eq. 
(5.26) respectively. These differences originate from the fact that we use n- 
dimensional regularization rather than n-dimcnsional reduction. The rcsTilt 
for the latter case can be obtained via a finite renormalization. It turns 
out that for n-dimensional reduction the rational numbers become equal 
to 7/4 for all three reactions which provides us with a strong check on our 
calculations. Finally we want to comment on the scale yU. In the computation 
of the radiative corrections we have assumed that the renormalization scale 
lir is equal to the mass factorization scale /x. If one wants to distinguish 
between both scales one has to substitute 




(5.27) 



in all finite expressions. 
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6 Differential distributions for the process 

In this section we will present the differential cross sections for Higgs-boson 
production in proton- proton collisions at the LHC. Here we study the de- 
pendence of the cross sections on input parameters like the QCD scale A, 
the renormalization/factorization scale fi and the dependence on the chosen 
set of parton densities. We also study which region in the gluon density is 
the most important. Furthermore we make a comparison with similar results 
presented in previous papers. Finally we also study at which values of the 
kinematical variables the soft-plus-virtual (S-l-V) gluons start to dominate 
the cross sections. This will give an indication about the validity of the S+V 
gluon approximation which is of importance when one wants to improve the 
perturbation series via resummation techniques leading to a better prediction 
for the differential cross sections. 

The hadronic cross section da is obtained from the partonic cross section 
daab as follows 

a J- a U a,b=q,g "'^l.min ''"1 "'^2,min 3^2 

X/«^(X2, //^) s'^^is, t, u, m^ n') .(6.1) 

In analogy to Eq. (2.12) the hadronic kinematical variables are defined by 

S^iPi + P2)\ T^iP,+p,)\ U^iP2+p,)\ (6.2) 

where Pi and P2 denote the momenta of hadrons Hi and H2 respectively (see 
Eq. (2.7)). In the case parton pi emerges from hadron Hi[Pi) and parton 
P2 emerges from hadron H2{P2) we can estabhsh the following relations 

pi ^ XiPi, P2 = X2 P2 , 

s — X1X2 S , t — Xi{T — m^) -|- , u — X2{U — m^) -\- , 
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When the matrix element behaves as I/S4, which occurs in the case of gluon 
bremsstrahlung or colhnear fermion pair emission, it is more convenient to 
choose the integration variable S4 instead of X2 in order to get better numer- 
ical stability. The cross section becomes 

^2 ^HiH2 .1 «4,^^, 



a 1 a U a,b=q,g -^a^i.min Xi Jo 



S4 — Xi{T — m?) — m? 



a t a u 



(6.4) 



The function X2{s4) can be derived from 

S4 — s + t + u — n? — X1X2 S + Xi{T — m^) + X2{U — m^) + , (6.5) 
and it reads 

xl{s,) = ^^-^^^T -ji^)- ^" ^ ^ ^^^s + T-m') + U. (6.6) 

xib + U — 

When parton pi emerges from hadron (-P2) and parton p2 emerges from 
Hi{Pi) one obtains the same expressions as in Eqs. (6.3)-(6.6) except that 
T and U are interchanged. This result has to be added to Eq. (6.4). When 
the partonic cross section is symmetric under t ^ u one can also use the 
representation in Eq. (6.4) without adding the result where T and U are in- 
terchanged provided one makes the replacement f^'^f^'^ — > f^^f^^ f^^f^^ ■ 
The expression in Eq. (6.4) simplifies when the partonic cross sections 
correspond to the Born reactions in Eq. (2.13) or the S+V contributions in 
Eqs. (5.24)-(5.26). In this case (Pa ~ ^(54) and the integral in Eq. (6.4) be- 
comes one-dimensional. For two-to-three body processes with no soft gluons 
or colhnear fermion pairs in the final state the integral is two dimensional as 
given in Eq. (6.4). If the processes are of the gluon bremsstrahlung type, 
or contain colhnear fermion pairs in the final state we have to split the two- 
to-three body partonic cross section in Eq. (4.40) into two parts i.e. a hard 
gluon part with S4 > A and a soft gluon part with S4 < A. The terms 
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containing the cut off parameter A cancel between the S+V cross sections 
originating from Eqs. (5.24)-(5.26) and the corresponding terms in the hard 
gluon hadronic cross section in Eq. (6.4), namely 

/' , ff' M^) ff'i^lM,,^) 

Jxi,min xi J A S4 — xi{T — m^) — m-' 



^2 ^HARD 

a t a u 



(6.7) 



In order to achieve the cancellation analytically we rewrite the expression 
above as follows 



/■I dxi 



■^4, max // S I TT r» TT o 



a t a u J A —xi[T — m^)—m^ 

j2 hard 

s4->o at a u 
•^a;i,min a;i JA —Xi{T — m^)—m'^ 

^2 ^HARD 

X lim s^^^, (s, t, S4, m^ /x^) . (6.8) 

s4_>o dtdu 

The expression between the square brackets is integrable in S4 so we can put 
A = in it. One can perform the integration over S4 analytically in the 
second part by using the simple expressions d'^a^^^^ in the limit S4 0. 
The ln*A//x^ terms, which arise from the integration over S4, then cancel 
those appearing in the S+V cross sections. 

In practice one is not interested in the Higgs boson differential cross 
sections depending on T and U but rather in the rapidity y and the transverse 
momentum px distributions. Neglecting the masses of the incoming hadrons 
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we have the following relations 

T — — y/S yjp^ + w? cosh y + ^Jp^ + w? sinh y , 

U — — y/S ^Jp\ + w? cosh y — ^fS \jp\ + sinh y , (6.9) 

so that the cross section becomes 

r/2 rT"iH2 ^2 H1H2 

5^-4— (5, y , m^) = 5^ (5, T, [/, m^) . (6.10) 
a pj. a y aT a U 

The kinematical boundaries are 



-b - r + <U < - 

from which one can derive 



m^-S<T<0, -S -T + m^ <U < ^ + , (6.11) 

T — 



IS 15 
< Pt < p|,i„ax , -T^ln— < 2/ < -In — 

2 2 



with p2 -G^±^!^_ 2 (g_;^2) 



or 



2 (S-m''y 2 

-|/max < 1/ < 2/max , < < ^-^ = p^^^^ , 

with t/ina^ = - In . , sg = .g,^ • 6.13) 

Since the cross section diverges for — we cannot perform the integral 
over this kinematical variable down to zero. However we can perform the 
integral over the rapidity and obtain the transverse momentum distribution 

d 0-^1^2 .j/^ax ^ (^2 ^HiH2 



dy- [S,pl,y,m'), (6.14) 

-3/max d pt d y 

with j/max given in Eq. (6.13). There is an alternative way to obtain the 
distribution above. This is shown in Appendix B. We checked that both 
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procedures lead to the same numerical result. In the case we plot the rapidity 
distribution we have to impose a cut pr,min on the transverse momentum 
integration i.e. 



d a 



H1H2 



d y 



iS,y,m') 



PT.max 
PT.min 



^2 cr"iH2 

dPT-j :j—{S,pl,y,m'^), (6.15) 

d Pt d y 



with pT.iimx given in Eq. (6.12). Actually the differential cross section dies 
off so fast that it is sufficient to put PT,ma.x = 8 x pT,mm- Finally we define 
what we mean by leading order (LO) and next-to-leading order (NLO). In 
LO the differential cross section is defined by 



rf2 a 



LO 



d Pt d y 



{S,PT,y,m'^) 



_d^a^ 
dpT dy 



{S,PT,y,m'^) , 



(6.16) 



and the gluon-gluon-Higgs coupling is given by G in Eq. (2.2) with C = 1 
(see Eq. (2.6)). We also adopt the leading logarithmic representation for the 
running coupling and the parton densities. The NLO corrected differential 
cross section reads 



^2 ^NLO 



a 



dpT dy 



{S,pT,y,m'^) = 



1 + 22 



dpT dy 




dpT dy 



{S,pl,y,m'^) 



{S,PT,y,m'^ 



(6.17) 



where the LO contribution in this formula is now multiplied with — 
1+22 as/ 4:7r. Furthermore the running coupling and parton densities are rep- 
resented in next-to-lcading order for which we have chosen the MS-scheme. 
In this way one obtains a result which is consistently corrected up to NLO. 
In our computations the number of light ffavours is taken to be = 5 which 
holds for the running coupling, the partonic cross sections and the number 
of quark fiavour densities. Further we have chosen for our plots the parton 
densities obtained from the sets MRST98 [32] CTEQ4 [33], GRV98 [34] and 
MRST99 [35] (sec Table 1). Notice that the GRV sets do not contain charm 
and bottom quark densities. For simplicity the factorization scale fi is set 
equal to the renormalization scale Hr- For our plots we take /x^ = + p^ 
unless mentioned otherwise. Here we want to emphasize that the magnitudes 
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MRST98 (LO, lo05a.dat) 


A^^ = 


: 130.5 MeV 




= 0.125 


MRST98 (NLO, ft08a.dat) 


^NLO 


= 220 MeV 


«f^°(M^) 


= 0.1175 


CTEQ4 (LO, cteq41.tbl) 




= 181 MeV 




= 0.132 


CTEQ4 (NLO, cteq4m.tbl) 


^NLO 


= 202 MeV 




= 0.116 


GRV98 (LO, grv981o.grid) 


A^o. 


= 131 MeV 




= 0.125 


GRV98 (NLO, grvnlm.grid) 


^NLO 


= 173 MeV 


a,^^°(Mz) 


= 0.114 


MRST99 (NLO, corOLdat) 


^NLO 


= 220 MeV 


a^^O(Mz) 


= 0.1175 



Table 1: Various parton density sets with the values for the QCD scale A 
and the running coupling ag- 

of the cross sections are extremely sensitive to the choice of the renormaliza- 
tion scale because the effective coupling constant G ~ «s(/Xr), which implies 
that da^'^ ~ and dcr^^^ ~ a^. However the slopes of the differential dis- 
tributions are less sensitive to the scale choice if they are only plotted over a 
limited range. 

For the computation of the Higgs-gluon-gluon effective coupling constant 
G given in Eq. (2.2) we choose the top-quark mass rrit = 173 A GeV/c^ and 
the Fermi constant Gp = 1.16639 GeV"^ = 4541.68 pb. In this paper we will 
only study Higgs boson production in proton-proton collisions at the center 
of mass energy y/S — 14 TeV characteristic of the LHC. Since the hadrons 
Hi and H2 are now identical the y differential cross sections are symmetric. 

The effect of the NLO corrections to Higgs-boson production were already 
studied earlier by the authors in [14]. In order to compare with their results 
we present LO and NLO differential cross sections in pt, integrated over y, 
for m = 120 GeV/c^ and //^ = + in Figs, la and lb respectively. 
The MRST98 parton densities [32] were used for these plots. We note that 
the NLO results from the q{q)g and qq channels are negative at small Pt so 
we have plotted their absolute values multiplied by 100. It is clear that the 
gg subprocess dominates but the g(g)g'-subprocess is also important. The 
corresponding results for the y distributions integrated over the pr region 
between j9T,min = 30 GeV/c and pT,max = 240 GeV/c are given in Figs. 2a 
and 2b respectively. The latter value was chosen because the cross section 
above 240 GeV/c is extremely small and can be neglected. Here the scale 
is fj,^ — m? + ]?r,min- Again we see that the g'g'-subprocess dominates. Using 
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the recent MRST99 set [35] we present the dependence of the NLO Pt- and 
l/-distributions on the Higgs mass in Figs. 3 and 4 respectively. 

Next we plot the scale dependence of the above differential cross sections 
in Figs. 5 and 6. For the pT-distributions we have chosen the scale factors 
^ = 2yU0) = A^o and = ^o/2 with /Xq = + p"^. In the case of the y- 
distributions we adopted Pq = +Pt mm where pT,mm = 30 GcV/c (sec Eq. 
(6.15)). Furthermore we have again adopted the MRST98 parton density set 
in [32] since it contains both LO and NLO versions. In the case of the px- 
distribution in Fig. 5a one observes a small reduction in the scale dependence 
while going from LO to NLO. This reduction becomes more visible when we 
plot the quantity 

N (ptA) = ^^j^^^^ (6.18) 
V IIqJ da{pT,lJ>o)/dpT 

in the range 0.1 < n/ Ho < 10 at fixed values of pr = 30, 70 and 100 
GeV/c. The upper set of curves at small /////o are for LO and the lower 
set are for NLO. Notice that the NLO plots at 70 and 100 arc extremely 
close to each other and it is hard to distinguish between them. Further one 
sees that the slopes of the LO curves are larger that the slopes of the NLO 
curves. This is an indication that there is better stability in NLO, which was 
expected. However there is no sign of a flattening or an optimum in either 
of these curves. This implies that one will have to calculate the differential 
cross sections in NNLO to find a better stability under scale variations. In 
the case of the ^/-distributions we adopted the scale /Xq = + PT,mm where 
PT,mm = 30 GeV/c (see Eq. (6.15)). As shown in Fig. 6 there is hardly any 
reduction in the scale dependence for the ^-distributions between the LO and 
NLO curves. 

Besides the dependence on the factorization and renormalization scales 
there are two other uncertainties which affect the predictive power of the 
theoretical cross sections. The first one concerns the rate of convergence of 
the perturbation series which is indicated by the X-factor defined by 

K = %^- (6.19) 



Another uncertainty is the dependence of the cross section on the specific 
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choice of parton densities, which can be expressed by the factors 

J CTEQ J GRV 

pCTEQ ^ " pGRV ^ " ^ (a nn\ 

The quantities in Eqs. (6.19) and (6.20) are plotted for the px distributions 
in Figs. 7a and 7b respectively where we have chosen the same parameters 
and scales (/x = /xq) as in Figs. 5a,b and 6. Depending on the parton density 
set the X-factors shown in Fig. 7a are pretty large and vary from 1.4 at 
Pt = 30 GeV/c to 1.7 at p^ = 150 GeV/c. Here the CTEQ4 parton densities 
lead to the smallest T^T-factor whereas the MRST98 set provides us with the 
largest one with the results from the GRV98 set in between. In Fig. 7b we 
show the dependences of the ratios defined in Eq. (6.20) as a function of pr- 
From this figure we infer that both the GRV98 and the CTEQ4 densities lead 
to larger cross sections than those computed from MRST98. The difference 
between the results obtained from the latter set with respect to the other 
ones is smaller in NLO than in LO. Furthermore In LO there is a small 
decrease in the R values as a function of pt which is in contrast to NLO 
where we observe an increase. In the case of the ^/-distributions in Fig. 8a 
again the CTEQ4 parton densities yield the smallest i^'-factor which does 
not vary much as a function of y. The latter also holds for the GRV98 set 
where however the fT-factor is larger than the one obtained from CTEQ4. 
The largest i^T-factor is obtained from the the MRST98 densities which show 
a much stronger dependence on the rapidity y than the ones obtained from 
the other sets. It becomes maximal at y = and decreases at larger absolute 
values of the rapidity. In Fig. 8b we show the ratio R as a function of the 
rapidity. Like in Fig. 7b the difference between the cross sections obtained 
from the MRST set and the two other sets becomes smaller while going from 
LO to NLO. At y = the discrepancy between the cross sections attains a 
maximum in the case of LO whereas it reaches a minimum for NLO. Notice 
that our results for the 7^-factors computed in NLO in Figs. 7a, 8a agree 
with those shown in the corresponding Figs. 2a, 2b in [14]. The same also 
holds for the NLO /^-factors in Figs. 7b, 8b. when compared with the same 
figures in [14]. The authors of [14] informed us that they did not adopt a 
fixed scale at //q = PT,mm + but they used a variable scale //q = Pt + 
when integrating over p^ to calculate the rapidity ^/-distributions. We have 
also run our programs with the latter scale, which yields different rapidity- 
distributions. The rapidity cross sections become smaller due to the smaller 
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running coupling constant but the i?-ratios hardly change w.r.t. the ones 
shown in Fig. 8b. Therefore the agreement between our results and those 
obtained in [14] are not spoiled. We also made a comparison between our 
results obtained for the differential cross section d a/d pt and those shown 
in Fig. la of [14]. To compare we have read off the central values of the bins 
from their Fig. la and, after changing bin sizes, have rcplottcd their values 
versus ours in Fig. 9 for the MRST98 set and the same input parameters. 
We only show the case where the scale is /j, — /Iq. Fig. 9 shows a small 
difference in LO and NLO between our cross sections (indicated in Fig. 9 by 
RSN(l)) and the those given by [14] (indicated by FGK). However the slopes 
are exactly the same. The difference might be due to a different choice of the 
effective Higgs-gluon-gluon coupling constant in Eq. (2.2). We have chosen 
rrit — 173.4 GeV/c^ whereas the authors in [14] took the limit — > oo (see 
Eq. (2.5)). If we take the value rrit — 10^ GeV/c^, which is rather close 
to an infinite top quark mass, our curves (indicated in Fig. 9 by RSN(2)) 
approach the ones given in Fig. la of [14]. For further checks we also made 
a comparison with the LO pT-distribution in Fig. 6 of [5] by reading off 
their values. Choosing the same input parameters and parton density set 
we completely obtain the same results as in [5]. Finally we checked several 
of the figures in [22]. In the latter reference only the NLO corrections to 
the 5'(?-subprocess have been completed. Using the same parton densities 
and parameters we found excellent agreement between our results for this 
subprocess and those obtained in [22] . 

We now analyse why the MRST98, GRV98 and CTEQ4 parton densities 
yield different results for the differential cross sections. This can be mainly 
attributed to the small x-behaviour of the gluon density because gluon-gluon 
fusion is the dominant production mechanism. To investigate the small x- 
behaviour we consider the product of the gluon flux with the corresponding 
partonic cross section 

F,,{xS,pl,fx') = %,{x,fx') ^{xS,pl,fx') (6.21) 

a Pt 

In this case we have removed the factor ag/^n from the LO (Born) cross 
section and the factor (Q;s/47r)^ from the NLO contribution to the par- 
tonic cross section. This was done to suppress the dependence on the strong 
coupling constant «s so that the differences between the several Fgg can be 
only attributed to the various parton density sets. In Fig. 10a we have 
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plotted the LO partonic cross section d cr^gg /d px versus log^Qa; for different 
values for pt- For x < 10^^ this cross section rises steeply whereas it flattens 
out for X > 10~^. The latter behaviour is changed when we look at the NLO 
partonic cross section d a^J /dpT plotted versus log^Q x in Fig. 10b. Here the 
rise at small x becomes even steeper whereas for x > 10^"^ the flatness shown 
by the Born cross section is replaced by a steep decrease when a; — >■ 1 except 
at high Pt- Notice that the NLO partonic cross section becomes negative 
(here for x > 8.0 x 10"^) due to mass factorization as explained below Eq. 
(5.15). 

To show the effect of the gluon flux wc have plotted Fgg{xS^p1,m^) as 
a function of Xog^^x in Fig. 11 for pt = 100 GeV/c using three parton 
density sets in table 1. From Fig. 11 we infer that the GRV98 set contains 
the steepest gluon density whereas the gluon densities from the other sets 
are about equal. From this observation one can understand the relative 
ordering of the plots in Figs. 7a and 7b. Another feature is that the hadronic 
cross section, which is represented by the integral of Fgg over x, receives 
its main support from the small x-region. Using the representation for the 
hadronic cross section d a/d pr as given in Eq. (B.7) we have computed this 
quantity where now all subprocesses are included. Further we have chosen 
two different integration regions. The first one is the full range Xmin < 
X < ajinax with Xinax = 1 (scc Eq. (B.8) and the second range is given 
by a;max = 5 x Xmin- In Fig. 12 we have shown the pT-distributions due 
to the two different integration regions. They do not differ by more than 
10% which shows that the whole differential cross section is dominated by 
the small x-region a; < 5 x Xmin. Hence apart from the different coupling 
constants the difference between the cross sections is due to the different 
gluon densities. Future HERA data will have to provide us with unique 
gluon densities before we can make more accurate predictions for the Higgs 
differential distributions. Note that we have chosen the axes in Fig. 12 to 
coincide with those in Fig. lb so that one can see the pt dependence of the 
i^T-factor for these two densities by overlaying the plots. 

After having studied the small x-region we now investigate the large x- 
region of the partonic cross sections where S+V gluons and coUinear quark 
anti-quark pairs dominate the radiative corrections. The S+V gluon part 
of the partonic cross section is obtained by omitting the hard contributions 
which are regular at S4 = so that we only keep the singular parts pre- 
sented in Eqs. (5.16)-(5.20). Furthermore we include the S-l-V partonic cross 
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sections in Eqs. (5.24)-(5.26). These two contributions constitute the S+V 
gluon approximation. To study its vaUdity we compute in NLO the ratio 



d 



^^^^ ~ J ^EXACT ' (6.22) 



for the Pt and y distributions which are given in Figs. 13 and 14 respectively. 
Here we use the MRST99 parton density set. One expects that the approx- 
imation becomes better at larger transverse momenta where pt — PT,max in 
Eq. (6.13) at the boundary of phase space, which leads io x — 1 using Eq. 
(B.8). However in Fig. 13 the highest value of p^, given by px = 150 GeV/c, 
is still very small with respect to PT,max ~ ^/S /2 = 7 x 10"^ GeV/c. There- 
fore it is rather fortuitous that the approximation works so well for px > 
100 GeV/c where one obtains it!^+^ < 1.2. The boundary of phase space is 
also approached when the Higgs mass increases. This can also be inferred 
from Eq. (B.8) where at fixed px, a^min ~^ 1 when S — pr- There- 

fore we have plotted R^^^ for various Higgs boson masses in Fig. 13. This 
figure reveals that the largest Higgs mass leads to the worst approximation 
contrary to our expectations which means that the kinematics is not in the 
large x-region. In Fig. 14 we plot the ratio R^'^^ for the ^-distributions for 
the same mass range as in Fig. 13 where we have chosen the cut pT.min = 30 
GeV/c in Eq. (6.15). Here the approximation is less good than for the pt 
distributions in Fig. 13 and like in the latter figure it becomes better when 
the mass of the Higgs gets smaller. At the lowest value chosen for the mass 
i.e. m — 120 GeV/c^ we observe that ~ 1.3. Furthermore we see 

no variation in with respect to the values of y. The approximation 

becomes better, see Fig. 15, if pr,min is chosen to be larger than 30 GeV/c 
so that it becomes closer to the values given in Fig. 13. Finally the fig- 
ures discussed above show that the S-l-V gluon approximation overestimates 
the exact NLO result. However this overestimate becomes smaller when the 
transverse momentum gets larger. In particular for p^- > 200 GeV/c the S+V 
approximation is good enough so that resummation techniques can be used 
to give a better estimate of Higgs boson production corrected up to all or- 
ders in perturbation theory. This statement will also hold when Higgs boson 
production is described according to the standard model approach where the 
boson is coupled via top-quark loops to the gluons without taking — > oo 
as we did above. Although the Pt distributions will change for the exact 
and S-l-V gluon approximation at large pr, the ratio R^^^ in Fig. 13 will be 
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less affected by the large top quark mass approach because of cancellations 
between the numerator and denominator in Eq. (6.22). Notice that the S+V 
gluon approach in the case of the exact cross section can be obtained from 
Eqs. (5.16)-(5.20) and (5.24)-(5.26) by replacing the approximate Born ma- 
trix elements M^jl^^. ^j, obtained from the effective Lagrangian in Eq. (2.1), 
by the exact expressions presented in [4], [5], [6]. 

Summarizing the above we have calculated the NLO corrections to the 
differential cross section for Higgs boson production in the large top quark 
mass approach which can be obtained from an effective Lagrangian. The 
calculation was carried out in standard n-dimensional regularization where 
the MS scheme was chosen for renormalization and mass factorization. We 
have presented some of the NLO results in the region where the effective La- 
grangian should be trustworthy. It turns out that the gluon-gluon subprocess 
dominates the hadronic cross section but the (anti-)quark-gluon subprocess 
is certainly not negligible. For the transverse momentum distributions there 
is a small reduction in scale dependence while going from the LO to the 
NLO cross section which is not observed for the rapidity distributions. Fur- 
ther there is still a large uncertainty in our predictions because the X-factor 
varies from approximately 1.4 to 1.7 depending on the parton density set. 
Also the dependence on the running coupling and the parton density set is 
appreciable. The latter is mainly due to the small x behaviour of the various 
gluon densities since both the partonic cross sections and the gluon densities 
increase very steeply at decreasing x. Finally we have shown that the S-l-V 
approximation is quite reasonable provided pT,min > 100 GeV/c in spite of 
the fact that x is still too small to belong to the large x-region. This means 
that this approximation can be used to rcsum the large corrections due to 
S-l-V gluons in order to obtain a better extimate of the all order corrected 
cross section. 
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Here we list the soft contributions to the partonic cross sections as defined 
in Eq. (4.42) which were not exphcitly given in section 4. For the various 

processes we obtain 
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j2 -soft 
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Appendix B 



There is an alternative way to obtain the transverse momentum distribution 
as presented in Eq. (6.14). First one computes the partonic cross section 
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The reason is that if one changes the variables t and u into S4 and p^ one 
has two possibilities 



t = 
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u 



or 
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which follows from the substitution 
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in Eq. (B.2). Hence one has to compute the sum dcr(ti, Ui) + da{t2, U2) which 
is equal to da{ti, ui) +da{ui, ti). Notice that when the partonic cross section 
is symmetric in t and u one can replace the sums above by 2 da{ti, Ui). This 
does not apply to the quark-gluon subprocess because here the cross section 
is asymmetric in t and u. The hadronic cross section is now obtained from 



(B.7) 

with 

r . = — V ^ -' 1 r = 1 fB 8) 

and denotes the partonic flux defined by 

= C dx, [' dx2 S{x - X, X2) f^'{x,,fx') f^%X2, . (B.9) 

JO ^0 

The lower boundary Xmin follows from the inequality in Eq. (6.13) where 
s = X S. We have checked that expressions Eqns. (6.14) and (B.7) lead to 
the same result. 
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Figure Captions 

Fig. la. The differential cross section d a/dpr integrated over the whole 
rapidity range (see Eq. (6.14)) with m = 120 GeV/c^ and fi^ = ni^+p^. 
The LO plots are presented for the subprocesses gg (long-dashed line), 
q{q)g (dot-dashed line) and 100 x {qq) (dotted line) using the parton 
density set MRST98(lo05a.dat). 

Fig. lb. Same as Fig. la in NLO except for 100 * abs(gg) (dotted line) and 
the additional subprocess 100 x abs(gg) (short-dashed line) using the 
parton density set MRST98(ft08a.dat). 

Fig. 2a. The differential cross section d ojdy for pT,min = 30 GeV/c and 
PT,ma^ = 240 GeV/c (see Eq. (6.15)) with m = 120 GeV/c^ and \j? = 
rn^ +PTmm- The LO plots are presented for the subprocesses gg (long- 
dashed line), q{q)g (dot-dashed hne) 100 x q{q) (dotted hne) using the 
parton density set MRST98 (lo05a.dat). 

Fig. 2b. Same as Fig. 2a in NLO with the additional subprocess 100 x qq 
(short-dashed line) using the parton density set MRST98(ft08a.dat). 

Fig. 3. The mass dependence of d a^^'^/dpT (Eq. (6.14)) using the set 
MRST99 with fj,^ ^ + p% for Higgs masses m = 120 GeV/c^ (solid 
line), m = 160 GeV/c^ (dashed hne) and m = 200 GeV/c^ (dot-dashed 
line) . 

Fig. 4. The mass dependence of d a^^^/dy for pT,min — 30 GeV/c (Eq. 
(6.15)) using the set MRST99 with /i^ = +Pr,min- The notation is 
the same as in Fig. 3. 

Fig. 5a. The scale dependence of d a^'^/dpr integrated over the whole ra- 
pidity range (see Eq. (6.14)) with m = 120 GeV/c^, /Xq = m^+p^ using 
the MRST98 parton density sets. The results are shown for ^ = 
(dashed line), = /iq (solid line) and = 2 x /^o (dot-dashed line). 
The upper three curves are the NLO results. 

Fig. 5b. The quantity N (pt , fJ. / fM)) (see Eq. (6.18)) plotted in the range 
0.1 < n/no < 10 at fixed values of pt with m = 120 GeV/c^ and 
— rn^ + p^ using the MRST98 parton density sets. The results are 
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shown for px = 30 GcV/c (solid line), Pt = 70 GcV/c (dashed line), 
Pt = 100 GeV/c (dot-dashed line). The upper three curves on the left 
hand side are the LO results whereas the lower three curves refer to 
NLO. 

Fig. 6. The scale dependences of d a^'^/dy and d a^^'^/dy integrated over 
the region pT-^min = 30 GeV/c andpT.max = 240 GeV/c (see Eq. (6.15)). 
The notation and parameters are as in Fig. 5. 

Fig. 7a. The X-factors in Eq. (6.19) for d a/dpr integrated over the whole 
rapidity range (see Eq. (6.14)) with m — 120 GeV/c^ and //^ = m^+p^ 
for the MRST98 sets (solid line), the GRV98 sets (dot-dashed hue), and 
the CTEQ98 sets, (dashed line). 

Fig. 7b. The ratios R in Eq. (6.20) for the same differential cross sections 
as in Fig. 7a, i?^^^ in LO (solid line), R^'^'^ in NLO (dot-dashed line), 
^CTEQ LO (dashed line) and i?CTEQ j^lO (dotted line). 

Fig. 8a The X-factors in Eq. (6.19) for d cr/dy integrated over the region 
PT,miii = 30 GeV/c and pT,max = 240 GeV/c (sec Eq. (6.15)) with 
m = 120 GeV/c^ and /x^ = +PT,min- The parton density sets and 
notation are as in Fig. 7a. 

Fig. 8b. The ratios R in Eq. (6.20) for the same differential cross sections 
as in Fig. 8a. The notation is the same as in Fig. 7b. 

Fig. 9. Comparisons of d a/dpx in LO and NLO obtained from our cal- 
culation (RSN) versus those in Fig. la in [14] (FGK), for rrit = 
173.4 GeV/c2 RSN(l) (sohd line), = 10^ GeV/c^ RSN(2) (long 
dashed hue) and rrit — oo FGK (dot-dashed line). The upper curves 
are for NLO and the lower ones for LO. 

Fig. 10a. The partonic differential cross section d a^g^ /dpT (Born) with m = 
120 GeV/c^ and /x^ ^ m'^ + p"^ for pr = 200 GeV/c (solid line), pr = 
150 GeV/c (long-dashed line), pr — 100 GeV/c (dot-dashed line), pr — 
75 GeV/c (short-dashed hue) and pt = 50 GeV/c (dotted line). 

Fig. 10b. Same as in Fig. 10a but now for the NLO differential cross section 
d af^/dpT. 
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Fig. 11. The flux multiplied by the partonic cross section represented by 
the quantity Fgg{xS,p\,rn?) in Eq. (6.21) with pt = 100 GeV/c, m = 
120 GeV/c^ and //^ = + plotted as a function of log^g^ 
MRST98 (solid hues), GRV98 (dot-dashed hues) and CTEQ4 (dashed 
lines). The upper curves are for NLO and the lower ones for LO. 

Fig. 12. The NLO hadronic cross section d a/dpx obtained from Eq. (B.7) 
(all subprocesses included) with m = 120 GeV/c^ and /x^ = + Pt- 
The parton density set is GRV98 with x^ax — 1 (solid line) and with 
Xmax = 5 Xmin (dashcd line). 

Fig. 13. The ratio R^^^ in Eq. (6.22) for the px distributions using the 

set MRST99 with p^ = + p\ j^j^ and various Higgs masses given 
by m = 120 GeV/c^ (solid line), m = 160 GeV/c^ (dashed line) and 
m = 200 GeV/c^ (dot-dashed line). 

Fig. 14. The ratio R^^^ in Eq. (6.22) for the y distributions using the 
MRST99 set with fj,"^ = rn? + Pr^min PT,min = 30 GeV/c and various 
Higgs masses given by m = 120 GeV/c^ (solid line), m — 160 GeV/c^ 
(dashed line) and m — 200 GeV/c^ (dot-dashed hne). 

Fig. 15. The ratio in Eq. (6.22) for the y distributions using the 

MRST99 set with //^ = + ^^^^ and m ^ 120 GeV/c^, for vari- 
ous values of pT,min, namely PT,min = 100 GeV/c (solid line), PT.min = 
200 GeV/c (dashed line), j9T,min = 250 GeV/c (dot-dashed hne) and 

Pt, min — 

300 GeV/c (dotted line). 
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